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In Aligned Natural Inflation, an alignment between different potential terms produces an inflaton 
excursion greater than the axion scales in the potential. We show that, starting from a general 
potential of two axions with two aligned potential terms, the effective theory for the resulting light 
direction is characterized by four parameters: an effective potential scale, an effective axion constant, 
and two extra parameters (related to ratios of the axion scales and the potential scales in the 2—held 
theory). For all choices of these extra parameters, the model can support inhation along valleys (in 
the 2—held space) that end in minima of the potential. This leads to a phenomenology similar to 
that of single held Natural Inhation. For a signihcant range of the extra two parameters, the model 
possesses also higher altitude inhationary trajectories passing through saddle points of the 2—held 
potential, and disconnected from any minimum. These plateaus end when the heavier direction 
becomes unstable, and therefore all of inhation takes place close to the saddle point, where - due 
to the higher altitude - the potential is hatter (smaller e parameter). As a consequence, a tensor- 
to-scalar ratio r = O (10““^ — 10“^) can be easily achieved in the allowed Us region, well within the 
latest Icr CMB contours. 


I. INTRODUCTION 

The Cosmic Microwave Background (CMB) data are in 
excellent agreement with the predictions of slow-roll in¬ 
flation [T]. Models that lead to sufficiently long inflation 
are very sensitive to Planck-scale physics. Even higher 
dimension Planck-Mass suppressed operators, such as 
4>^ y/Mi, □ induce corrections of O (1) or larger to the 
slow roll parameter 77, and therefore shorten the duration 
of inflation to a few e-folds. This UV-sensitivity is present 
for any inflationary potential, but it particularly affects 
models that generate a relatively large gravity wave sig¬ 
nal (tensor-to-scalar ratio r > 0.01), as in such models 
the inflaton vacuum expectation value (vev) changes by 
an amount greater than the Planck scale during inflation 
0 - 

A well known and simple symmetry that can forbid 
such operators is a shift symmetry, namely the invariance 
of the action under the transformation </> —>■ ^-1-constant. 
A field possessing this symmetry, at least at an approx¬ 
imate level, is denoted as an axion. In order to have a 
nontrivial <()—dependent potential, the symmetry cannot 
be exact, but it can be broken at a technically naturally 
small level: assume that only one sector of the theory 


^ In this expression, 0 denotes the inflaton, V its potential, and 
Mp is the (reduced) Planck scale M-p = where Gw 

is Newton’s contant. 


breaks the shift symmetry, leading to a nontrivial poten¬ 
tial for (/i; if this potential is small, quantum corrections 
to it will also be small, as they must be proportional to 
it (loops that do not contain such term do not break the 
shift symmetry, and therefore do not generate a poten¬ 
tial for the axion). The first work that recognized the 
merits of employing a shift symmetry in inflation is [3], 
where the resulting class of models was named Natural 
Inflation. 

A relevant scale in axion inflation is the axion scale 
/. This is the scale that determines the least irrelevant 
shift-symmetric coupling, often the dimension five cou¬ 
pling ^FF to gauge fields. In presence of this operator, 
gauge instantons break nonperturbatively the continuous 
shift symmetry down to the discrete cj) —)■ (j)+2Trf. There¬ 
fore, / controls the periodicity of the axion potential, in 
absence of any explicit breaking of the shift-symmetry. 
Refs. [Hll] studied the simplest potential with such prop¬ 
erty 


U = A^ 


1 — cos 
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As conventionally done in the literature, in this work we 
use the term Natural Inflation exclusively for Q, but 
we credit mm for realizing the relevance of the shift- 
symmetry protection, and not simply of the potential 

0 - 

Models of axion inflation are under control for 
m^, H < f, where m,p denotes the axion mass, and H 
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the Hubble rate during inflation. This is due to the fact 
that the theory of the axion is obtained by integrating out 
modes that are heavier than /. Ideally, one would also 
want to have f < Mp (see [S] for a partial list of relevant 
literature on this issue). One reason for this is that one 
expects quantum gravity to break the shift symmetry, as 
any global symmetry, at the Planck scale. This argument 
does not apply if the shift symmetry originates from a lo¬ 
cal symmetry, as typically in string theory. However, all 
known controlled strong theory constructions have sub- 
Planckian axion scales. The / < Mp request contrasts 
with the fact that the model Q is consistent with the 
CMB results only for f > 7 Mp [T] . Due to the more 
speculative nature of these arguments, the conservative 
approach of relaxing the / < Mp requirement can be jus¬ 
tifiable. Nonetheless, a large number of studies have been 
appeared in the recent literature, with the goal of pro¬ 
viding a UV completion of natural inflation, formulated 
in terms of axions with a sub-Planckian scale. 

Mechanisms developed in these works include: iden¬ 
tifying the axion with the extra component of a gauge 
field in a 5d theory compactified on a circle EHSl ; us¬ 
ing two axions with a nearly aligned potential m, using 

= O (lO^ — 10^) axions El in the mechanism of as¬ 
sisted inflation |12j : generating a monodromy through 
an explicit breaking of the shift symmetry [UIE5] ; slow¬ 
ing down the inflaton from particle production, HZl- 
IdO] (as in the original model of warm inflation EH); 
coupling the axions to higher forms [32H84] : consider¬ 
ing a non-standard mixing of the axion with gravity 
[35| : modifying the axion evolution with a coupling to 
a non-abelian vector field with nonvanishing vev |3tij , 
considering higher derivative terms [421143] ; considering 
multiple sinusoidal functions |lil446j : realizing an align¬ 
ment from kinetic mixing |47H49j ; explicitly breaking the 
shift-symmetry via (fMFF operators [50l[5T]; employing 
Coleman-Weinberg contributions to the axion potential 
[521153] . Alternatively, models with non-compact axions 
have also recently been considered [Ml [SS] . 

Very recently, a discussion has started in the literature 
of whether also models with sub-Planckian axions that 
produce an effective super-Planckian excursion (as the 
model of uni that we study here) are also affected by 
large quantum gravity corrections [56H61j . An argument 
in this discussion is that fulfilling the Weak Gravity Con¬ 
jecture (WGC) [52] in such models implies the existence of 
gravitational instantons that produce unsuppressed cor¬ 
rection to the potential. The conditions under which such 
corrections are relevant, and their magnitude, are still 
under debate. In particular, ref. [501 discussed condi¬ 
tions under which there can be gravitational instantons 


^ A related mode, that can be understood as the limit of m in a 
certain parameter range, has been formulated in |37| . Both the 
original models of m and m are ruled out by either stability 
considerations or the CMB data [381440| . See also ref. m for a 
recent extension of m to two axions. 


that fulfill the WGC, without introducing unsuppressed 
higher harmonics to the potential. We do not address 
these aspects in the present work, but we rather study 
the phenomenology of the aligned model of m taken at 
face value. Our goal is to show that the model possesses 
inflationary trajectories that have not been emerged in 
the many studies of the model, and that have a particu¬ 
larly interesting phenomenology (as, contrary to the pre¬ 
dictions obtained from Q , they are compatible with the 
CMB data at Icr). Although our computations are re¬ 
stricted to cni, we hope that inflationary solutions with 
analogous properties can also be present in other mod¬ 
els of multi-field axion inflation, and that our result can 
motivate the search for such solutions. 

Prior to the present work, the alignment mechanism of 
[TU] was studied in [551IM] , and extended to the N > 2 
axions case in [55H59] . The works [751185] studied the em¬ 
bedding of the alignment mechanism in string theory. 
Despite this extensve series of works, a detailed study of 
the phenomenological predictions of Aligned Natural In¬ 
flation [T5| is still missing. This is the goal of the present 
work. 

The most general potential for the two axions intro¬ 
duced in m is characterized by two potential terms and 
four axion scales. An alignment between these scales re¬ 
sults in a hierarchy between the mass scales of the two 
axions. We denote the light and heavy state by (j) and 
ip, respectively. Despite the simplicity of the model, ob¬ 
taining the phenomenology for all possible choices of the 
parameters requires extensive computations. We improve 
over existing discussions of the mechanism by computing 
the inflationary trajectories ip {<p) in the model, and the 
effective 1—field potential along such trajectories. We 
find that, in the limit of strong alignment, this effective 
potential is characterized by four parameters: an effec¬ 
tive axion scale an overall potential scale A^, and two 
additional parameters, rf and ta, related to ratios of the 
axion scales and the potential scales of the starting po¬ 
tential. For any fixed rf and r\, the model has the same 
freedom as Natural Inflation: the value of controls 
the overall flatness of the potential (and so the amount 
of inflation along an inflationary trajectory), while the 
scale A is fixed by the normalization of the scalar per¬ 
turbations. Therefore, we learn that, to study the gen¬ 
eral phenomenology of the Aligned Natural Infaltion, we 
“only” need to vary over the combinations r/ and 
in addition to what is done for the single field model of 
Natural Inflation. 

An inflationary trajectory pj {(p) needs to have (i) suf¬ 
ficient flatness along the light direction, and (ii) stability 

T r 

along the heavy direction, > 0. From our analysis 
we learn that the trajectories of Aligned Natural Infla¬ 
tion behave in qualitatively different ways according to 


^ Other studies of axion inflation in string theory can be found in 
[8111821 f84U92| . For a generic study of inflationary models based 
on alignment, see m- 
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the relative value oi r\ vs rf. For all choices of these two 
parameters, the model admits inflationary trajectories 
connected to minima of the potential: the heavy direc¬ 
tion is stable along these trajectories, and inflation ends 
close to the minimum, when the (f> direction is no longer 
flat. These trajectories lead to a phenomenology anal¬ 
ogous to that of Natural Inflation. As for this model, 
the predicted values of and r lie on a region which 
is outside the latest 1 a CMB contours [T] . For a sig¬ 
nificant range of the extra parameters (specifically, for 

J'Ae 3 in the case of r/> 1, and for rAC 4,4 in 

the case of r/ < 1) we find the presence of a new class of 
trajectories that has not emerged in previous discussions 
of the mechanism. This new class of trajectories occurs 
near a saddle point of the full potential, and inflation ter¬ 
minates not because a minimum has been reached, but 
because the heavy direction becomes unstable. Thanks 
to this feature, the flatness of these trajectories along the 
light direction is decoupled from the end of inflation, and 
inflation occurs on a higher altitude plateau, where the 
potential is flatter. This leads to a much smaller value 
of r with respect to Natural Inflation, well within the Icr 
CMB bounds. 

Another characteristic of models of axion inflation is 
that the shift symmetry highly constrains the coupling 
of the inflaton to matter. One expects that the domi¬ 
nant among such interactions is the j(j)FF coupling to 
gauge fields, where c is a dimensionless coefficient that 
in the spirit of effective field theories is expected to be 
of O (1). As first shown in [95], a coupling % > 


/ ~ IQie GoV 

of the inflaton to any light gauge field results in overpro¬ 
duction of gauge quanta during inflation. In the limiting 
case j = O ( iQie^sev ) a production can have inter¬ 
esting phenomenological consequences, that we review in 
Section |V| While this production is negligible in Natural 
Inflation (since f > 7 Mp is required), a large enough cou¬ 
pling can be easily present in Aligned Natural Inflation. 
We update the discussion of this effect already presented 
in the literature by computing the precise value of the 
coupling necessary to obtain a sizable effect along some 
of the inflationary trajectories obtained in this work. 


The paper is organized as follows. In Section jn] we 
review the model introduced in m, and we obtain the 
effective 1—field potential in the limit of strong align¬ 
ment. In Section jlTlj we discuss more in details the light 
field trajectories in the 2—field potential. In Section jlV] 
we study the resulting phenomenology in the {rig — r} 
plane. In Section jV] we discuss possible additional signa¬ 
tures of Aligned Natural Inflation that can emerge from 
the couplings of the axions to light gauge fields during 


inflation. In Section VI we present our conclusions. In 
Appendix we present some analytic properties of the 


This behavior is reminiscent of hybrid inflation [Mi- 


inflationary trajectories. In Appendixj^we briefly review 
the computation of cosmological perturbations in 2-field 
inflationary models, that we use to obtain some of the 
results of Section |iy| In Appendix [Cj we warn against an 
approximation that does not produce the correct 1—field 
description of the model. 


II. NATURAL ALIGNED INFLATION AND 
THE STRONG ALIGNMENT LIMIT 


The general 2—field {9, p) model of Natural Aligned 
Inflation is characterized by four axion constants and two 
potential scales [TO] 


V = K 


4 

1 


1 — COS 






1 — cos 




( 2 ) 

This potential has one exact flat direction, which we de¬ 
note by 0, for fi/gi = f 2 /g 2 - If this relation is only 
approximately valid, that direction becomes flatter than 

A2 j ^2 

the naive expectation, namely ^ y , ^. Equiv¬ 
alently, this can be seen as having an axion (j) with an 
effective scale » fi, gi- This is the aligned mecha¬ 
nism of m- 

To study this model, we find covenient to express the 
two potential scales as 


A 


4 _ 
1 — 


1 

1 + T’A 


A^ 


AN 


r-A 

7 + rA 


A^ 


(3) 


so that A is an overall potential scale, and ta is the ratio 
between the two scales in ([^. We also find convenient to 
express the ratio between the axion scales as 


4. = ^ = ^3 

fi 1 +a ' /2 1 - a ’ 


Finally, we redefine 


fi=rff, /2 = — . 


(5) 


In this way, we reformulate the model (§ in terms of 
two parameters of mass dimension one (/ and A) and 
four dimensionless parameters: 


A4 1 + 


cos 


ta 


1 

ff \f 


1 -I- a 

rg l) _ 


ta 

- cos 

1 + J'A 





1 - « p\ 
rg f)_ 


( 6 ) 


The potential is symmetric under = 

|’"A ^ ^ ^ ^ ly } (iiiferchange of the 

two terms) and under S 2 = {9 ^ —9 , p ^ —p} (the 
potential is even). We stress that ^ is equivalent to 
(i). and that no approximation has yet been made. 

In this work we obtain phenomenological results for 
Aligned Natural Inflation using the exact 2—field model 
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^ and the effective 1 —field model that we derive be¬ 
low. The excellent agreement between the results from 
the 2 —field and the 1 —field model confirms the accuracy 
of the approximation that we perform in the reminder of 
this Section. Specifically, we are interested in the model 
in the limit of strong alignment. In our parametrization, 
this means |a| <C 1. In the following, we study the the¬ 
ory to the leading nontrivial order in a. We start by 
performing the rotation 


= R 


R^ 1 — R-i-i — 


■0 


I-rirA 


^1 + rj (l + r 2 )"/'l + ^/^A 


1-rir A 


Ri 2 = -R21 = - a ■ • \ " .( 7 ) 

^l + r2 (l + r-2)'/"l + ^/^A 


The rotation matrix is orthogonal up to O (a^) correc¬ 
tions, and therefore the kinetic term in the new basis 
reads 


[dOf {dpf 

2 2 


(M" 

2 


M! 

2 


[l + 0(a2)] ,(8) 


and in the following we disregard the O (a^) corrections. 

The transformation Q rotates (up to O (a^) correc¬ 
tions) the fields 0, p in tne light </) and heavy ^ eigenstate 
of the potential ([^ expanded at quadratic order around 
the minimum located at the origin. Along the inflation¬ 
ary trajectories, ip is always the heavy field (up to O (a) 
corrections), and the trajectories are obtained by inte¬ 
grating Ip out. Integrating ip out leads to a trajectory 
Ip = Ips ( 0 ), which has the property ips = O ( 0 ^) next to 
the origin, but that is ips ^ Q in general. It is crucial to 
compute the precise trajectory ips ( 0 ) in order to obtain 
the correct 1 —field effective description of the model. 

Under the rotation 0 . the potential becomes 


_U 

A4 


Cl 


C2 


C 3 


C 4 


1 - 


cos (ci a 0 -I- C2 Ip) rA cos (03 a 0 -I- C4 ip) 


1 + CA 

2rjrA 




1 + rirA 


1 + CA 

[1-1-0 (a)] , 


l + r2 


f^f Tg 


[1 + 0 (a)] 
2rf 


/y^l + c2 (l + rj ta) 


T [1 + C)(«)] 


frg 


[1 + 0 (a)] 


( 9 ) 


and in the following we disregard the O (a) corrections to 
the coefficients q. This expression for the potential man¬ 
ifestly shows that 0 is the light direction in the limit of 
large alignment, since it enters multiplied by a <C 1 and 
by an O factor in the argument of the two cosines, 

while Ip enters times an O ^ factor without any a sup¬ 
pression. n Therefore, apart from special points where 
= 0 , the curvature of the potential in the ip direction 
is much greater than the curvature in the 0 direction. 

To visualize and study the potential, it is convenient 
to perform a different rescaling of the 2 —fields: 


2 a 0_0 7 + ^9 0 _ Ip 

-7 = W ’ ^ ^-7 = W ’ 

7/0 ^g j /b 


l + r2 


This allows us to draw a contour plot of the potential 
V ^0, 0^ in which the curvature appears to be compara¬ 
ble in both directions, and the trajectories can be better 
visualized. However, one should keep in mind that, in 
the basis { 0 , 0 } of the canonically normalized fields the 
curvature is much greater in the ip direction. As a con¬ 
sequence, starting from a generic initial condition, the 
system first evolves in the ip direction, with practically 
no change in 0 , until it reaches a minimum ^ = 0 . 
Then, the system continues to evolve along a trajectory 
that minimizes V with respect to ip. 

In terms of the rescaled fields (10), the action of the 
model is 



rfijj (. 11 ) 


This formulation of the theory is characterized by 5 pa¬ 
rameters, and it is equivalent to the original theory up 
to 0(a) corrections. The potential in (11) is still in¬ 
variant under the symmetries 51^2 defined above, which 

now read Si = |r-A ^ 7 y, 0 ^-- 0 ,r/^ 7 )^| and 

'S '2 = {0 -t -0,0 -t -0|- 

In the lagrangian the fact that 0 is a much lighter 
field than ip is encoded in the kinetic term, given that 

f2 

^ = O (a^) +; 1. Therefore, in the classical field evolu- 

tion, the kinetic energy associated to 0 can be neglected 
with respect to that associated to 0 (it would provide 


® In Appendix]^ we show explicitly through some examples that 
setting -0 = 0 everywhere provides in general an incorrect CMB 
phenomenology. 


We are implicitly assuming that, apart from a, the other dimen¬ 
sionless parameters in ([ 6 |l are of order one. 
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a correction to the kinetic term of the same order as 
those that we have neglected in eq. (§)• Disregarding 
the kinetic term of ip in 0 amounts in treating ip as a 
lagrange multiplier. We can then integrate ip out: 


dV 

—^ = 0 solved by ip = ips 
dip 


( 12 ) 


This solution provides the light field trajectory in the 
2—field space; this trajectory is stable provided that 
> 0. When this is the case, the trajectory can 
support inflation provided that is large enough. In¬ 
serting the solution (121 back into 0 we obtain the one 
dimensional effective theory 


£=|(9 




Vi 


= V2 


tps 


(13) 

The 1—field effective theory (131 has 4 parameters: the 


overall potential scale A, the effective axion scale fcj,, 
and the two dimensionless parameters Vf and re¬ 
lated to ratios between the axions scales and the po¬ 
tential scales of the original theory, respectively. Due 
to the symmetries Si and S 2 , the 1—field effective po¬ 
tential is separately invariant under cp ^ —(p and under 
^ , ta—T herefore we can restrict our dis¬ 
cussion to r/ > 1. Specifically, in the main text we focus 
our attention to r/ > 1. The limiting case r/ = 1, which 
is one of the very few cases for which we could solve (12) 
analytically, is discussed in Appendix [X] 

We conclude this Section by further discussing the 
2—field potential (11). It is easy to verify that all the crit¬ 


ical points of the potential take place when both terms 


in (111 are at an extremum. If we denote by A and B 
the argument of the two cosines in 0, this happens 
when A = Ti n, and B = tt m, with integer n and m. 
The 2—field potential has a minimum when both m and 
n are even, a maximum when they are both odd, and 
a saddle point when they are one even and one odd. 
Moreover, due to the symmetry of the potential under 

A ^ A + 2tt , B —)■ B + 27r, we can divide the ip'^ 
plane in an infinite number of equivalent domains. 

In Figurej^we show one such domain, delimited by the 
lines ^ = 0, —27r and B = 0, 2tt. The domain has the 

origin, with "(/il coordinates 


O: {0,0}, 


(14) 


in one of its corners (the other three corners are three 
equivalent minima of the potential). The central part of 
the domain is occupied by the maximum 


M : TT < — 



(15) 


Finally, each of the four sides contains a saddle point, 
with equivalent saddle points on opposite sides. We de¬ 
note by Sa and Sn, respectively, the saddle poin t on the 



0 5 10 


FIG. 1. Contour plot of the potential as a function of the 
rescaled fields <p (horizontal direction) and ip (vertical direc¬ 
tion). Due to the different rescaling for the 2—fields, the 
potential exhibits comparable curvatnre in both directions. 
However, for |a| <C 1, the field ip is significantly heavier than 
(p, and the evolution proceeds along trajectories where ^ = 0 


^ = 0 and on the B = 0 line. These points are marked 
in Figure and they have coordinates 


Sa = tt 


1 1 

r/’ l + rjrAj 


Sb = TT 



rf 1 
i + rjrA j 
(16) 


Without loss of generality, we can restrict the initial 
conditions for the fields to be in this domain. Moreover 

- again due to the symmetry properties of the potential 

- we can restrict the initial conditions to be along the 
valley that ends on O, or (when it exists) along a valley 
that starts from one of the two saddle points Sa^b- Such 
valleys are shown in Figure and are studied in Section 
|III[ Any other valley in the potential can be mapped to 
one of these valleys. 


III. VALLEYS AND CRESTS OF THE 2-FIELD 
POTENTIAL 
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rf = \5\r^= 0.1975 r/ = 1.5; = 0.4144 



FIG. 2. Contour plot of the potential centered in the domain at the right of the origin. The a;—axis is 0 = while the i/—axis 


is '0 = ^, see eq. |l0| ). The green (magenta) curves are stable valleys (unstable crests) within the domain that pass through 
the origin, the maximum, and / or the saddle points Sa, Sb- The plots are for r/ = 1.5 and for different values of va in the 
characteristic regions: ta < \ (top-left); \ < r a < \ (top-right); ta = \ (bottom-left); rA> \ (bottom-right). 

r f r f r f r f r. 


We study the valleys and crests of the potential 0; 
these are lines in the plane for which = 0. 

The stable valleys minimize the potential in the '0 direc¬ 
tion > 0), while the unstable crests maximize it 

( < 0). As discussed at the end of the previous Sec¬ 

tion, we restrict our attention to rj > 1, and to valleys 
and crests that are (i) inside the domain we have singled 
out, and (ii) connected either to the origin or to one of 
the saddle points Sa,b- The equipotential lines of Figure 
show how the valleys and crests extend to other do¬ 
mains, or beyond the maximum Ai. However, due to the 
symmetries of the potential, we do not need to consider 
these extension in the present discussion. 

In general, we could not solve the condition ^ = 0 
analytically, so that we do not have general analytic ex¬ 


pressions for the valleys and crests tps ■ However, it 

is straightforward to obtain such lines numerically. We 
summarize the result of our investigation in Figure 
where valleys and crests are shown as green and magenta 
lines, respectively. The results summarized in the figure 
agree with the analytic results that we could obtain in 
neighborhoods of the extrema of the potential (the ori¬ 
gin, the maximum, and the saddle points), and for the 
special case ta = ^. This analytic results are reported 
in Appendix and they confirm the validity of what 
we present here. The valleys and crests have a qualita¬ 
tively different behavior according to the relative value 
of ta vs ry. We discuss these different behaviors starting 
from small ta, and discussing how the lines change as ta 
increases. 


For ta < ^ (top-left panel of Figure 


2), a valley ex- 
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r/ = 1.5; ta = 0.4144; rg = 1; a= 0.01; fi=53 Mp 


FIG. 3. Contour plot of |ll| , for parameters in the ^ < 

region, together with valleys (green) and crests (magenta). 

The two red curves are two distinct inflationary trajectories 
in this model. They are obtained from a numerical evolution 
of the exact model §. Both evolutions shown contain 60 e- 
folds of inflation plus a brief transient moment after inflation 
in which the system reaches the minimum. 


tends from the minimum O to the saddle point Sa, and a 
crest extends from the maximum M to the saddle point 
Sb- In the interval < t’a < ^ (top-right panel) we 

still find a first curve 'ips extending from the min¬ 
imum to Sa, and a second curve extending from the 
maximum to Sb- However both curves are partially a 
stable valley and partially an unstable crest. For ta = ^ 

(bottom-left panel), the two curves join each other at a 
point. Finally, for > r'j (bottom-right panel), a stable 
valley extends from the origin to Sb, while an unstable 
crest extends from the maximum to 5 a. 

At ta = ^, a reconnection between the two tps 

curves take place. Specifically, for < \ (top two 

panels of Figure^, a curve connects the origin with 5 a, 
and a curve connects the maximum with Sb- Instead, 
for ta > ^ (bottom-right panel of Figure 2), the origin 

is connected with Sb and the maximum with 5 a. A 
“swap” between two halves of the two curves takes place 
at TA = A-, which is the only value of ca for which the 

two curves touch each other at a point. 

Provided the effective axion scale is sufficiently 
large, the valleys are inflationary attractors. The val¬ 
leys shown in Figure are obtained from the potential 
(111, where O (a) subdominant terms in the exact poten¬ 
tial (§ are disregarded. We performed several numerical 
evolutions in the exact potential for a = and 


we verified that the valleys of Figure ([^ describe the 
inflationary evolution with very good accuracy. Figure 

shows two examples of this. We choose parameters 

in the -^ < r\ < \ interval, so that two distinct and 
’■/ ”/ 

inequivalent N valleys exist, one connected to the mini¬ 
mum at the origin, and one disconnected from it. The 
red curves in the figure show two inflationary trajectories 
obtained from the exact potential §. These are two dif¬ 
ferent inflationary evolutions for the same model and for 
the same choice of parameters, but with different initial 
conditions. As we discuss below, they lead to a different 
phenomenology. 

The trajectory on the inflationary valley disconnected 
from the origin is particularly interesting. Inflation ends 
not because the potential reaches a minimum, but be¬ 
cause the inflationary valley (^- > 0) becomes an un¬ 
stable crest (^- < 0), so that inflation is terminated by 
an instability in the heavy direction. We see in the fig¬ 
ure that, after leaving the inflationary valley, the system 
reaches a new = 0 line; the heavy field '0 performs 
damped oscillations about this valley, and eventually the 
system reaches the minimum on the top of the figure by 
evolving in the light (j) direction. This stage of the evo¬ 
lution is very fast. The evolution shown in the figure 
is characterized by iV = 60 e-folds of inflation along the 
valley connected to Sb', the following phase, from the mo¬ 
ment the system leaves the valley to when it first reaches 
the minimum, lasts for ~ 2.9 e-folds. During this sec¬ 
ond stage, the equation of state oscillates with average 
Wave — —0.16. This phase should be understood as the 
beginning of the post-inflationary reheating. 


IV. {us - r} PHENOMENOLOGY 


In this Section we study the CMB phenomenology of 
Aligned Natural Inflation in the {ug — r} plane. As dis¬ 
cussed in the previous Section, we find two classes of 
inflationary trajectories in this model: those along val¬ 
leys connected to a minimum, and those along valleys 
disconnected from any minimum. In the first case, in¬ 
flation ends as the fields approach the minimum of the 
potential; in the second case inflation terminates at the 
end of the valley, due to an instability in the heavy tp 
direction. This second class of solutions exist only for ta 


the 


1 1 
”/J 


interval 


0 


^ In the sense that the valleys cannot be mapped into each other 
by a symmetry of the potential, so they give a different phe¬ 
nomenology. 

® To be precise, such an evolution can also take place for r\ slightly 
greater than so that evolutions where inflation ends due to 

instability in the ip direction are possible for ^ < "Ta < ^ + e, 

with e small. The reason for this is that, for ta = one finds 













r 

0.14 ^ 



FIG. 4. Potential along the two valleys shown in Figure]^ 


Significantly different phenomenological results are ob¬ 
tained in these two cases, and, for this reason, we discuss 
them in two separate subsections. We anticipate that the 
evolutions along valleys disconnected from minima lead 
to a significantly smaller value of the tensor-to-scalar ra¬ 
tio r than the other trajectories, and of what is generi- 
cally found in axion inflation. The reason for this is that 
the potential is flatter ( = smaller e parameter) closer to 
a saddle point than to a minimum (see Figure]^. We 
stress that the valleys disconnected from the minima are 
missed if one rotates the potential at the origin into a 
heavy tp and light (p direction at the origin, and he/she 
then makes the approximation that all of inflation takes 
place along the ip = 0 line. 

We performed several inflationary evolutions for the 
model, with different values of the parameters, and we 
verified that the CMB phenomenology obtained from the 
1—field effective potential (13) is in very good agreement 
with that obtained from the exact model (|^. In the plots 
shown below we have obtained the fields evolution from 
the exact model (§, and we have then computed the 
values of Us and r using the slow-roll relations (B19). 


A. Inflation on valleys connected to minima 

For definiteness, we computed the CMB phenomenol¬ 
ogy for the case in which the axion scales are all compa- 


= 0 at the precise point where the valley connects with the 
crest shown in the figure, so that the fields do not bend along the 
valley, but escape from it (ending inflation), and again reach the 
minimum shown on the top of the figure. This behavior rapidly 
disappears as r\ increases slightly above since > 0 all 

along the valley in this case. For the parameters used in Figure 
jsj we numerically found that e ~ 0.004 (while ^ ~ 0.4444). 



FIG. 5. Predictions of Aligned Natural Inflation (with in¬ 
flation along a trajectory connected to a minimum) in the 
{ris — f} plane, confronted with the la and 2cr Planck con¬ 
tour lines (specifically, we choose to plot the more conservative 
red contour lines of Figure 12 of ID)- We fixed r/ = 1.5 and 
varied r\ as follows: the dashed lines, from bottom to top, 
are for ta = 0.01,0.07,0.1,0.19; the solid lines, from top to 
bottom, are for ta = 0.3, 0.4144, 0.5,1, 3. The lowest theoret¬ 
ical curve, drawn as a dotted line, is for Natural Inflation. All 
the theoretical curves are done for V = 60 e-folds of inflation. 


rable to each other, r/ = 1.5 and = 1. We also fixed 
the alignment parameter to a = and the number of 
e-folds of inflation to TV = 60. We then varied ta to cover 
all the different cases shown in Figure The results of 
our evaluations are shown in Figure Each line shown 
in the Figure corresponds to a given value of ta and to 
varying f^. 

All curves have a common end on the top-right part, 
corresponding to large f^. Figure]^ shows the potential 
along a connected valley as a function of ^ As 

grows, a given value of (p produces a longer inflationary 
expansion. Therefore, if N is fixed to 60, one needs to 
decrease the value of -p as increases. At large one 


is probing only the part of V 



that is closest to the 


origin, where a quadratic approximation of the potential 
suffices. One then recovers the values of massive chaotic 
inflation, Ug — 0.967 and r ~ 0.133. 

As a comparison, the black dotted line in the figure 
refers to Natural Inflation [3]. All the four dashed lines 
in the figure sample the model in the ta < 4 region. 

We see that the smallest value of ta shown gives results 
very close to Natural Inflation. This is explained by the 
analytic computation presented in Appendix where 
we show that Aligned Natural Inflation reproduces Nat¬ 
ural Inflation in the limit of very small r\. We then see 
that increasing r\ in this interval leads to progressively 
greater values of r. This is particularly true in the left 
portion of the curves shown in the figure. As we men¬ 
tioned above, the top-right part of the curves is obtained 
at relatively large while decreases as one moves to¬ 
wards the left of the curves. As decreases, ( p ne eds to 
start closer and closer to the saddle point. Eq. (A61 gives 
the analytic form of the effective I—field potential close 
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to this point. We see from this relation that increasing 
TA (in the ta < ^ regime that we are considering here) 

indeed results in a less flat potential, and so in a greater 
value of r, in agreement with the dashed curves of the 
figure. 

The two top solid theoretical curves shown in Figure 
sample the \ < rf^ < \ regime. We recall that in this 

figure we only consider inflation in the valley connected to 
the minimum. We see from the top-right panel of Figure 
[^that this valley becomes a crest before reaching a saddle 
point. For this reason, cannot be taken arbitrarily 
small, and, as a consequence, the corresponding curve 
in the {n^, r} plane only extends for a finite interval. 
The length of this interval decreases as we approach the 
r\ = ^ value. This is due to the fact that also the 

fraction of the stable portion of the curve from Sa to the 
minimum decreases as ta —> 

'■/ 

Finally, the three bottom solid theoretical curves 
shown in Figure ^ sample the ta > r'j region. At rel¬ 
atively small ftj, ^eft part of the curves), th e ev olution 
starts close to the saddle point Sb- Eq. (A8) shows 


that the potential becomes progressively flatter as in¬ 
creases in this regime. This explains the behavior of the 
three bottom solid curves in the figure. We see that the 
curve corresponding to the highest value of r\ shown also 
approaches that of Natural Inflation. This is explained 
by the analytic computation presented in Appendix 
where we show that Aligned Natural Inflation reproduces 
Natural Inflation in the limit of very large r\. 

We also performed computations for other values of ry, 
and we obtained results qualitatively similar to those of 
Figure For brevity, we do not show them here. 


B. Inflation on valleys disconnected from minima 

As we mentioned, in the \ < ta < ^ interval, infla- 

tion can occur on valleys that are connected to the saddle 
point Sb and that are disconnected from any minimum. 
Inflation ends because the heavy direction becomes un¬ 
stable. The tensor-to-scalar ratio in these trajectories is 
significantly smaller than the one found for the valleys 
connected to a minimum. This is visible in the two Fig¬ 
ures]^ and where, respectively, the two cases rf = 1.5 
and Tf = 8 are studied (such values do not have any 
particular importance, and they have been chosen just 
as a representative case of comparable axion scales, or 
somewhat hierarchical axion scales). 

It is possible to reproduce analytically the results 
shown in these two figures with good accuracy. Most 
of inflation occurs close to the saddle point, where the 
1—field effective potential reads 

0 < ^ < ^0 ■ (17) 



r 



FIG. 6. Predictions of Aligned Natural Inflation (with in¬ 
flation along a trajectory disconnected from a minimum) in 
the {hs — r} plane, confronted with the Icr and 2(j Planck 
contour lines. The theoretical lines have been obtained 
for Tf = 1.5, and for va, from bottom to top, equal to 
0.25,0.33,0.38,0.41,0.43. All the theoretical curves are done 
for A = 60 e-folds of inflation. 

r 



FIG. 7. Predictions of Aligned Natural Inflation (with in¬ 
flation along a trajectory disconnected from a minimum) in 
the {ug — r} plane, confronted with the Icr and 2cr Planck 
contour lines. The theoretical solid lines have been ob¬ 
tained for rf — 8, and for ta, from bottom to top, equal 
to 0.005,0.009,0.012,0.014,0.015. The dotted line visible in 
the top-left corner is the theoretical prediction of Natural In¬ 
flation. All the theoretical curves are done for A = 60 e-folds 
of inflation. 


In this equation, we have defined (cf. eq. (A8)), 


‘ = Trr//^ -(j) , / = ^ W - 

r/ y TA 


(18) 


while (j)Q denotes the value of (j) at which inflation ends 
(due to the fact that the orthogonal direction ip becomes 
tachyonic at that point). This is a top-hat potential with 
maximum on the saddle point Sb (of the 2—field poten¬ 
tial) located at (/> = 0. The slow roll parameters on this 
potential are 




rj ~ 


2M2 


<^ « / • 


(19) 


Integrating the slow roll relation ^ ~ one ob¬ 

tains the value of at A e-folds before the end of infla- 
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tion: 


(j)M f" 


( 20 ) 


and finally, from the slow roll relations r ~ 16e , — 1 ~ 

2r} — Qe ~ 2ry, one obtains the functional relation 


r ~ 8 (1 - Tis) 


PQ^_N(l-n,) 


P 


( 21 ) 


Each theoretical curve in Figures and is obtained 
for fixed rf and r\, and for varying f^. The value of 

^ q 2 xr 

(po can be obtained from imposing = 0 along the 

inflationary trajectory, where the potential 0 is used. 

We see that any fixed value of r/ and ta, results in a 
12 

given value of ^. Therefore, each curve is characterized 
by the same value r of so that (21) gives the pre¬ 


/' 


diction for the shape of the theoreticaTcurves, up to an 
overall rescaling. A comparison with the exact numerical 
results shown in Figures and shows that (21) - with 
00 obtained from a numerical evolution with the exact 
2—field potential - is accurate up to a few percent level 
discrepancy for ta — ^, while the discrepancy increases 

to ^ 10 — 50% level for ta — 

We can have a fully analytic relation by estim ating 
the value 0o at which inflation vanishes. Eq. (A7) is an 


analytic expression for the inflationary trajectory in the 
neighborhood of 5^. We insert this expression in 

expand to quadratic order in i50, and set this expression 
to zero. This provides a simple estimate for the value of 
0 at which the heavy direction becomes tachyonic: 


00 

/ 




- 1 


v^-v/1 - ta Pf 


( 22 ) 


We stress that this value cannot be extremely accurate, 
as it is obtained by expanding the potential close to 
Sb^ and not close to the where the instability occurs. 
Nonetheless, it provides a reasonable estimate. By com¬ 
paring with the exact numerical evolutions used for Fig¬ 
ures [^and[^ we find that the estimate ( [2^ is accurate up 
to ~ 10% discrepancy for ta — while the discrepancy 

increases to ^ 20 — 30% level as ta —t Combining 


(21) and (22), one obtains 


- l) 

1 - TA 


(1 _ „^) e-A^(i-"») . (23) 


Despite not perfectly accurate, eq. (23) can be used to 


understand the behavior of r with the parameters of the 


9 


This is confirmed with very good accuracy by the numerical evo¬ 
lutions with the exact potential ||m. 


potential. For instance, it shows that, at any fixed rf, 
the value of r grows with growing r\, in agreement with 
what is seen in Figures and It also shows that r in 
general grows with rf (t he d erivative with respect to rf 
of the right hand side of (23), evaluated for the r\ > \, 


is positive). This is also confirmed by the comparison 
between Figure and Figure 


V. COUPLING TO VECTOR FIELDS 

The couplings of an axion to “matter” fields are highly 
constrained by the shift symmetry. The leading operators 
that give the coupling to gauge and fermion fields are, 
respectively, ^ (p F F and 0757'^V’: where / is the 

axion scale, and where ca and are model dependent 
dimensionless coefficients that, in the spirit of an effective 
field theory, are naturally expected to be 0(1). 

The coupling to gauge fields typically controls the per¬ 
turbative decay of an axion inflaton after inflation, as 
the decay into fermions is helicity suppressed [S]. More¬ 
over, this coupling also leads to interesting gauge field 
production already during inflation. C3 Specifically, one 
gauge field helicity becomes tachyonic in presence of this 
operator m, and the corresponding gauge quanta are 
non-perturbatively produced (this effect vanishes both 
in the UV regime and IR regime, so the production is 
“self-regulated”, and does not diverge). The amount of 
produced quanta is exponentially sensitive to the pa¬ 
rameter [57] ^ , in the C > 1 regime, while the 

production is negligible at ^ < 1 (for definiteness, we 
assume ^ to be positive in this discussion. Changing 
sign of ^ simply changes the gauge helicity state that 
is produced, and so it does not change the phenomeno¬ 
logical constraints). These gauge quanta inverse decay 
to produce inflaton perturbations, that are highly non- 
gaussian |95j and that have a blue spectrum that can 
lead to too many primordial black holes [lOOj . They can 
also produce gravity waves |96j that violate parity unii, 
and that also have a blue spectrum, growing to a level 
that can be tested already at the next generation inter¬ 
ferometer experiments |98l I102L I103j . All these effects 
constrain ^ < 0(1). The estimate of |100j suggests that 
the limit from primordial black holes, that forces ^ < 1.5, 
is the strongest one, although the precise limit is sensitive 
to the evolution of the scalar perturbations in a highly 
nonlinear regime, which may not be fully under compu¬ 
tational control nnsj. The other limits are more robust, 


These considerations are valid for any gauge field whose quanta 
are light during inflation. Although explicit computations in the 
literature assume U(l) fields, the results can be easily extended 
to the non-abelian case. See [5] and m for the details of the 
computations that we discuss in this section. 

In this discussion we consider the production of gauge fields from 
the inflaton during inflation. See m for a discussion and m 
for a more detailed study of the production at reheating. 
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and constrain ^ < 2 — 2.5 0.0 Using the slow roll 
relation \^\ ~ H Mp, we can rewrite 


c = 



(24) 


and so we see that the regime ^ > 0(1) is typically 
reached when — < 10“^Mj, 1951 . 

Such effects may help in phenomenologically distin¬ 
guishing the single field Natural Inflation from the 
Aligned Inflation mechanism. Indeed, as we now discuss, 
they are negligible in the former case, while they may be 
detectable in the latter case. To see this, consider the 
operators 

AUNatural — ~ ~ ^J ^ ^ ^ 

AUAligned = ~ ~ ^ 

where / in the second line is the axion scale introduced 
in ([^. The parameters c,C 0 ,Cp are model-dependent di¬ 
mensionless coefficients (for instance, such a coupling can 
arise from integrating out heavy fermions that couple 
both to the axion and to the gauge field), which may 
naturally be expected to be of order one. In the aligned 
inflation case, we use Q to rotate to the base of light 
{(j)) and heavy field (tp) during inflation, and we arrive to 
the coupling 


AUaligned — ^ 


o^fcp 


c^tp 

u 


FF 


(26) 


where the scales and have been defined in (10), 
while, again, the parameters and are model- 
dependent dimensionless coefficients which may naturally 
be expected to be of order one. ^ 


As we have discussed after eq. ( [10| ), the rescaled fields 
(p = -^ and '0 = ^ exhibit a similar excursion during 
inflation. However, the original field (p moves much mor e 
than Ip. This is encoded in the ^ factor entering in (26): 
due to the greater excursion, the light field (p is cou^d 
to the gauge field much more strongly than ip. So, the 
coupling of the infiaton to the gauge field in Natural Infla¬ 
tion and Aligned Natural Inflation is given, respectively, 
by 


^natural — C ^ 



^aligned 



We therefore see that, in Natural Inflation, a value c > 
O (lO^ — 10^) is necessary to produce a visible effect, in 



FIG. 8. The quantity is a dimensionless coefficient con¬ 
trolling the coupling of the infiaton in Natural Aligned Infla¬ 
tion to a light gauge field (defined in eq. (26l). The curves 
show that value of this coefficient leading to observable gauge 
quanta production during inflation, for the special choice of 
parameters r/ = 1.5 and ta = 0.4, and for an inflationary 
trajectory along a valley that is connected to / disconnected 
from a minimum of the potential. 


contrast to the natural c < 0(1) expectation. Also in 
Aligned Natural Inflation, a value ^ ^ O (lO^ — 10^) 
is required to have a visible effect. However, due to the 
- 1 enhancement, this can be possible even if < 

0 ( 1 ). 

In Figure 8 we show the required value of to pro¬ 
duce ^ = 1 m Aligned Natural Inflation, for the param¬ 
eter choice Tf = 1.5, va = 0.43. We observe that a 
stronger coupling is required for trajectories along val¬ 
leys disconnected from minima of the potential, since 
the infiaton potential is flatter along those valleys (see 
the discussion at the beginning of Section IV). 

To summarize, in Natural Infiaton the request of a suf¬ 
ficiently fiat potential, namely / > Mp, imposes a strong 
limit on the coupling strength j x 0(1) of the infiaton 
to gauge fields, which precludes an interesting produc¬ 
tion of gauge quanta during inflation. On the contrary, 
if the aligned mechanism is employed, one can start from 
significantly smaller axion decay constants, and hence a 
significantly stronger coupling to gauge fields, and still 


obtain a fiat infiaton potential, /(/> = O > Mp. The 

fact that none of the effects associated to the gauge field 
production has been yet observed translates into limits 
of the coupling of the axions to any light gauge field, and 
provides useful constrains for model building [98) . 


The quantity ^ typically increases during inflation. The limits 
quoted here refer to the value assumed by ^ when the large CMB 

scales left the horizon. _ 

These coefficients are related to those entering in l |25[ ) by = 

——and C.^ = , up to O (o) corrections. 


Figure|^has been obtained from the exact value of p} in numerical 
evolutions of the exact 2—held model § , with a = 0.01 and 
Tg = 1. However, the same result can also be obtained from the 
effective 1—field model ( |13| |, where the parameters a and Tg do 
not enter (one only requires that a <C 1). 
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VI. CONCLUSIONS 


The alignment mechanism of [TO] has received con¬ 
siderable attention in the recent past. Despite of this, 
a detailed study of the inflationary trajectories in this 
model, and of the associated phenomenology, was still 
missing. The main goal of this work was to perform such 
a study. We were particularly interested in understand¬ 
ing whether the model can be phenomenologically dis¬ 
tinguished from single field Natural Inflation. The pre¬ 
diction from Natural Inflation are outside the current la 
CMB contours, and the model will be fully explored once 
the sensitivity to r ~ 0.03 is reached m, see for in¬ 
stance Figure (^0 A number of current and near future 
CMB experiments have claimed r = O (l0“^) sensitivity, 
and the current discussion of Stage 4 experiments aims 
to reach r = O (l0“^) sensitivity |106) . It is therefore 
of particular interest to study whether Aligned Natural 
Inflation can provide values of — r below the dotted 
line shown in Figure as this region cannot be covered 
by single field Natural Inflation. 

The most interesting result emerged from our analysis 
is the existence of inflationary trajectories near saddle 
points of the potential, and disconnected from minima of 
the potential (in the sense that these trajectories end far 
away from minima, because the heavy held perpendicular 
to them becomes unstable). As studied in Section IV 
these trajectories lead to much smaller values of r than 
those obtained in Natural Inhation. The Ug — r prediction 
for the evolutions along these trajectories is given by eq. 
(23), which is a fully analytic expression in terms of the 


parameters in the model and of the number of e-folds. 
As we discussed, this relation is only an approximate one, 
but it is reasonably accurate, and it allows to understand 
the dependence of the phenomenological predictions on 
the parameters of the model. 

A second possible way to distinguish Aligned Natu¬ 
ral Inhation from single held Natural Inhation originates 
from the fact that a stronger coupling of the inhaton to 
other helds should be expected in the former, given that 
the starting axion scales are sub-Planckian. A moder¬ 
ate level of alignment, a = O(l0“^), is associated to 
O (lO^® GeV) axion scales, where the coupling of the ax- 
ions to any massless or light (to < H) gauge held can 
lead to signihcant gauge held production during inhation, 
with several phenomenological consequences [5]. We re¬ 
viewed this in Section FV! 

Finally, we hope that the conclusions learnt from the 
present study can be of relevance also for other multi- 


In our figures, we fixed the number of e-folds of visible inflation 
to Nqmb = 60. Such a value is obtained for a quick reheating 
after Natural Inflation. Slower reheating (under the conventional 
assumption that the equation of state at reheating is smaller than 
that of radiation ESS], 'i'^reh ^ §) implies that the large scale 
CMB perturbations were produced closer to the end of inflation 
(VcMB < 60), when the potential was less flat. This gives a 
larger value of r than the one shown in Figure [sl 


field models of axion inflation. Recently, there has been 
considerable interest of inflation in the landscape of string 
theory, and on the relevance of critical points in multi- 
field inflationary potentials, see for instance |107H112] . 
We have shown that the presence of critical points in 
one of the 2—held extensions of Natural Inhation can 
strongly impact the phenomenological predictions of that 
model, and this may be the case also for other multi- 
held extensions that we have mentioned in this work. 
It may be interesting to perform a systematic study of 
trajectories, and explore the relevance of saddle points, 
also for such models. 
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Appendix A: Analytic properties of the valleys and 
crests 


In this Appendix, we denote as “curve” a line ips 


along which = 0. Namely 




sm — 


1+rjr^ 


, 2 
H- +rf TA sm 

) 


1 + rjrA 


+ Tf Ipsj = 0 . 

(Al) 


These curves are stable valleys if > 0, in which 

case they can support inhation provided that the effec¬ 
tive axion scale is sufficiently large. They are instead 
unstable crests if < 0. 


Inserting the solution of (Al) into the 2—held 

potential V 2 provides the effective 1—held potentia l, as 
formally described in eq. (13). In general, eq. (Al) 


can be solved numerically. Here we discuss the cases for 
which we could obtain an analytic solution. 


1. Critical points 


The condition (Al) is obviously satished at the critical 
points O, A4, Sa, Sb listed in Section |n] We can there¬ 
fore solve it analytically for small displacements about 
these points. Next to the origin we obtain 


V's - - 


rjvA 


(r^TA + l) 


(A2) 


By computing and imposing (|A2|), we can ver¬ 


ify that (A2) is a stable valley (which is obvious, given 


that the origin is a minimum of ¥ 2 ). We observe that. 
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as the valley reaches the origin, 'ips vanishes faster than 
(j). This is due to the fact that, as we discussed after per¬ 
forming the rotation the 2—fields (j) and tp are mass 
eigenvectors near the origin, so that the light field valley 
must proc eed along ^ = 0 at this point. We also see 
from (A2) that, on the immediate right of the origin, the 


valley bends to positive ib for ta < -^ and to negative ip 
for rA > ^. This can be seen, respectively, in the top 
two panels of Figure and in the bottom right panel. 

Being of O , the field ip does not contribute to the 
quadratic potential near the origin, where we have 


stable near Sb for ta > 


TyT - - 

Expanding along (|A7|), we find that this curve is 




which is the case for the 


top-right panel and the two bottom panels of Figure 
Along the trajectory, and close to Sb, the potential reads 


bl^neartSB — 


1 


1 

1 - V 


r? TA 


4: rjrA — 1 


5(p^ 


(AS) 


2 . 


1"A = 


1 


lA,nearO — o 

2 (1 


r^TA 


■rA) ( 


1 + rirA 


(A3) 


Near the maximum, namely for ^ = + Sp 


and ips = TT 
by 


Tf TA-rf 


Sip, the condition (All is solved 


Sip — 




6 


rj rj - 1 


(r^TA + l) 


sp^ 


(A4) 


This curve is unstable near P4. To the left of the max¬ 
imum, this curve bends towards negative p for ta < ^ 

’■/ 

and towards positive p for > \- This is the oppo- 

site of the curve connected to the origin. As Figure 
shows, when one curve connects the minimum with the 
saddle point Sa at positive p, the other curve connects 
the maximum with the saddle point 5 b at negative p, 
and viceversa. 

Near the saddle point 5 a, namely for 0 = fj + Sp and 


V's = i+rirA A condition (A1) is solved by 


2 ri ta 

^ 8^2 ■ 


(A5) 


Expanding along (A5), we find that this curve is 
stable near 5 a for r\ < -V, which is the case shown in 

the top-left panel of Figure Along the curve, and close 
to 5a, the potential reads 


El 


l.neartSA — 


2rA 
1 T^a 


1 

1 - T 


4 1 — rA 


Sp'^ 


(A6) 


Finally, near the saddle point Sb, namely for p = 
TTTf + Sp and ps = — TT-ir-h Sp, the condition (Al) 

v l-\-T 

is solved by 


2 ri rA 

^ 8^2 ■ 


(A7) 


As discussed in the main text, at the interesting value 
of ta = a reconnection between the two ps (pp^ 

curves takes place. For ta = ^, the field p enters in 


the two arguments of the two sine functions in (Al) mul¬ 
tiplied by the same coefficient . y^e could then solve 


(Al) analytically, and hnd that the condition is sat¬ 
isfied by two straight lines: one joining the minimum O 
with the maximum Ai, and one joining the two saddle 
points 5 a and Sb- The two lines intersect at the point 
C, of coordinates 


C-- ^ j — 
2 V/ 


r/, 0 


(A9) 


One can show analytically that only parts of the two 
curves are stable (^^ > 0): specifically, the segment 
between C and the minimum O, and the segment be¬ 
tween Sb and C. The other parts of the trajectories are 
unstable. 

Along the stable C — O segment, the potential reads 


El = 1 — cos 



(AlO) 


and the argument of the cosine ranges from ^ to 0 as one 
moves from C to O along the trajectory. 

Along the stable Sb — C segment, the potential reads 


|.2 _ 1 

El = 1 -I- - 4 , -- sin 

r2 + l 



(All) 


and the argument of the sine ranges from | to 0 as one 
moves from 5 b to C along the curve. 


3. rf = l 


We recall that the potential (11) is invariant under the 


simultaneous exchanges rj ^ ^ and ca ——. For this 
reason, we do not need to discuss the case r/ < 1. In the 
main text we studied the case Vf > 1. Here we discuss 
the rf = l case. The main difference with the discussion 
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of the main text is that now the regime ^ < ta < ^ is 

'’/ '■/ 

absent. 

For rf = 1, the field ip enters with the same coefficient 
in the two potential terms. Thanks to this, we can write 


1^2 = 


di) 

where 


A cos 'tp + B sin ■)/> 

1 + ?’A 

AUj smjp - B Uj cosp) Q 2 y^ 


1 + T'A 


dpP- 


= 1 - ^^2 , 
(A12) 


Aicp] = cos ^ —h ta cos 

1 + r-A 

D / A \ 

B { (p ] = sm ^-TA sm 


1 + J'A ’ 


1 + J'A 


1 +?’A 


(A13) 


For TA = 1, one finds A = 2 cos and B = 
0. Therefore, the condition ^ = 0 leads to the two 

o-ip 

curves '0s = 0 (a straight line connecting the minimum 
O to the maximum Ai), and 0 = tt (a straight line con¬ 
necting the two saddle points). Along the first curve, the 

potential reads Fl = 1 — cos and so the trajectory 
is stable in its first half close to the minimum). Along 
the second curve, the potential is constant, Fi = 1, and 

o2 T/- 

yrAr = 0. This can be understood as the limit rf —>■ 1 
of the case ta = ^ shown in the bottom-left panel of 
Figure 

Let us now discuss the ta ^ 1 case. We write the 
solutions of = 0 as the two curves satisfying 


sin '0s = ± 



/A2 


B^ 



(A14) 

The curve with the -I- sign connects the minimum O with 
the saddle point Sa for ta < 1, and with the saddle point 
Sb for TA > l.[^ In both cases the variable 0 ranges from 
0 to TT, and the effective potential along the trajectory 
evaluates to 


F = l- 


-I- -I- 2 ta cos 0 

1 + r-A 


(A15) 


Since the pote ntial i s always < 1, the curve is a stable 
valley (see eq. (A12)). 


16 ■\Yg recall that the model is symmetric under ry —t ^ and ta -A 
In the current case, ry = 1, the symmetry identifies ta 
with We see that the saddle points interchange role in this 
identification. 


The curve with the — sign instead connects the maxi¬ 
mum At to the saddle point Sb for ta < 1, and with the 
saddle point Sa for ta > 1. This curve in an unstable 
crest. The behavior and the stability of the two curves 
agree with those discussed in the main text for rj > 1 
(cf. the top-left panel of Figure for ta < 1, and the 
bottom-right panel of Figurej^for ta > 1), showing that 
the rf ^ 1 limit is continuous. 


4. Very small or large ta 


These are the regimes in which one term in the poten¬ 
tial ([^ is strongly dominant over the other term. The 
condition = 0 sets this term to a negligible value. 
Then one is left with the other term, evaluated along 
the trajectory set by the first term. This reproduces the 
potential of Natural Inflation. 

To be specific, the very small (respectively, large) ta 
regime requires ta < (respectively, ta > ^). 

We start the discussion with the very small ta regime. 
In this limit, the left hand side of eq. ( |Al[ ) can vanish 
only if the the argument of the first sine is at most of 
O 1- This implies 

4 


and, as a c onsequence, the argument of the second sine 
in eq. (A1) is ry 0 -|- O ■ Using this information 

in the potential (111, we learn that 

A4 


F = 


1 —cos^r/0^ -|-O (r^r^) . (A17) 


1 + rA L 

Namely, up to negligible corrections. 




1 — cos ^rf 0^ 


1 1 

, rA ^ ~2 ^ “4 

i n 


(A18) 


A completely analogous computation shows that 


F~ At 



1 

’’A > “2 ■ 


1 

04 


(A19) 


Appendix B: Cosmological perturbations in the 2 
field model 


We start by computing the cosmological scalar pertur¬ 
bations for a generic model of two canonically normalized 
real scalar fields, characterized by the action 

(Bl) 


























15 


We work in conformal time r and in spatially flat gauge: 

(t) [— (1 + 2(/)) dr^ + 2diBdx^dT + Sijdx^dx^] . 

(B2) 

Moreover, we decompose the scalar fields as a 
T— dependent background component plus perturbations 


‘fi = (jii (r) + (5(^i(T,f) = (j)i (r) + 


d^k 


^ Qi (k) 

^ik-x \ / 


(27r 


, 3/2 


a(r) 
(B3) 

We insert (B2) and (B3| into the starting action, and 
obtain the quadratic order action in the perturbations 
Si- We then integrate out the non-dynamical fields B 
and (j). □ Doing so, one obtains 

s, = Jdrd^kc , ^ ^ [\qt - fc'iftr - q:, 


2 

m2 = 

“ M2 


3 _ ^'k^'k 


a'2 2M2 j 


+ 


M2 O' 


; (0'D, + ) 


— av 


2M2 


2a 


l2-\ 


where V^i = and prime denotes differentiation with 
respect to r. We note that the variables Qi are the canon¬ 
ical variables of the problem, and that their conjugate 
momenta are 


n = ^ = Qt' 
* “ dQ[ 


(B5) 


We further decompose the Fourier modes Qi into mode 
functions times annihilation/creation operators: 

Qi (jtj = Vij (k) aj (jtj + V^j (fc) a] (^-kj ■ (B6) 

The mode functions satisfy the classical equations of mo¬ 
tion 


V" + n^V = 0 , = /c2 1 + M 


(B7) 


where 1 is the identity operator, and M^ the squared 
mass matrix defined in (B4). To quantize the system. 


we impose standard Equal Time Commutation Relations 
among the fields and their conjugate variables, and the 
standard algebra on the annihilation/creation operators: 


[Q^ (t, x ), Ilj {t, y)] = iSij (x-y) , 


(^) ’ 


= 6ij (jt - 


(B8) 


The fields 0 oc Sqqq and B oc Sgoi enter in S 2 without time 
derivatives, and so they do not add any extra dynamical degree 
of freedom. It is customary (see for instance \m) to integrate 
them out (i) by evaluating the constraint equations = 

0 , (ii) by solving them expressing the nondynamical variables (p 
and B in terms of the dynamical variables Qi and Q 2 -, and (iii) 
by inserting these expressions back into the action 82 - In this 
way 82 becomes an action in terms of the dynamical fields only. 


To ensure these relations, we impose 

ppt' = -i , 

VV^ = real , V' — real, 


(B9) 


at some given initial time tq. The first condition in (B9) 


ensures that both relations (B8) can be imposed at tq. 


Enforcing all three of (B9) at tq ensures that these con¬ 


ditions continue to hold at later times (as can be imme¬ 
diately verified taking their time derivatives and using 
(B7)). This ensures that the first relation in (B8) holds 


at all times. 

We impose the initial conditions when a mode is deeply 
inside the horizon. Specifically, the physical momen¬ 


tum 


a(To) 


should dominate the dispersion relation of 


(B4) 


the mode. The entries in the mass matrix (B4| are of 
O where mi _2 are the two mass eigenval¬ 
ues of M. Therefore, we ensure that H, m, . Ne- 

glecting slow roll corrections, a = during inflation, 
and so we impose 


^ ^Hro) = k^ 


1 + 0 


1717 


fc2 Tq ’ k'^H'^TQ yj 


(BIO) 

Using this approximated expression for D, we see that 
the equations of motion (B7) admit the early time solu¬ 


tion Djj,early - ■ So, we impose. 


(t-o) = 


\/2k 


^IJ , 


77', (to) = -i\ 7 ; k 


Hj 5 


(Bll) 


which satisfy all three of (B9|. Namely, in the deep sub¬ 


horizon regime the mass matrix M is irrelevant, and the 
two eigenmodes of the system behave as decoupled mass¬ 
less modes in Minkowski space-time. The solution just 
written is the standard adiabatic vacuum solution for 
these decoupled modes. 

We are interested in observables that are linear in the 
two scalar fields, 0{t,x) = Ci{t) S(j)i{t, x), where the 
(real) coefficients Ci can depend on background quantities 
(and, hence, on time). We assume that we have a set of 
such observables 0^'^\ ..., 0^"^! so that 


0 ^“^ {t,x) = Car (t) S(j)i {t, X) , 


(B12) 


and we compute the two point correlation function 
among pairs of them 


Cap {t, X, y) = (t, x) 0^^'> {t, y) 

f°° dk sin (kr) 


kr 


Pap {t,k) , r = \x-y\, 


(B13) 


where we have used the fact that (thanks to statistical 
isotropy) the mode functions V do not depend on the 
direction of A:, and where we have introduced the power 
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FIG. 9. Solid lines: Auto- and cross- power spectra for the 
adiabatic and entropy scalar perturbations in Aligned Natural 
Inflation, obtained from a numerical evolution with the model 
Q, for the parameters ta = 0.3, r/ = 1.5, = 1, = 

5%/2Mp, comoving momentum of the 

mode that exits the horizon 60 e-folds before the end of infla¬ 
tion. Inflation starts next to the saddle point Sb, and ends 
when the trajectory becomes unstable in the heavy direction. 
Dashed line: Slow roll result (B181 for Pniz in the single field 
approximation (assuming 5 = 0). 


spectra 


Pai 3 (r, k) = [c (r) V (r, k) (r, k) (t)] , 

(B14) 

(we note that this expression is symmetric and real). 

In shorts, for any set of observables (B12), we compute 


the self- and cross-power spectra by imposing the initial 
conditions (Bll), by evolving (B7|, and by evaluating 
dBlil). 


In this work, we have two scalar fields (pi = 0 and (p 2 = 
p, and we are interested in the power in the adiabatic TZ 
and entropy S modes, that, in spatially flat gauge, read 
(see |114j for a review) 


H f cos 9 sin 9 

|0/| \^—sin0 COS0 


5 ( 1)1 

S(j)2 


(B15) 


where we have defined 


COS0 = , sin0 = , \$'\ 

w\ w\ 


^ \N)? + 5)? ■ (B16) 


The solid lines of Figure show the auto- and cross¬ 
spectra of the adiabatic and entropy modes for an illus¬ 
trative choice of parameters, with inflation taking place 
along the high-altitude plateau connected to the saddle 
point Sb- The adiabatic mode corresponds to pertur¬ 
bations along the light inflationary direction, while the 
entropy mode corresponds to the much heavier orthogo¬ 
nal direction (the alignment parameter is a = j^). As 


expected, the entropy mode is significantly smaller than 
the adiabatic one, and we can simply disregard it, with 
very good accuracy. Starting f rom the quadratic action 
(B4), performing the rotation (B15l, and setting 5 = 0, 


we obtain the theory for the single TZ perturbation, in 
the limit in which S can be disregarded. As expected, 
we explicitly verified that, in this limit, the quadratic 
action for TZ reproduces that of a single field with a po¬ 
tential equal to the potential of the 2—field model along 
the inflationary trajectory. Specifically, setting 5 = 0, 
we obtain 


V _a\(j)'\ 

TZ= - , z= —— 
z H 


S 2 = - f drcPk 




■ (B17) 


which is the standard action ma for the Mukhanov- 
Sasaki variable [ml [nil in the single field case. 

One may then use the standard slow roll relations for 
the perturbations. 


Pn = 






- 1 = 2?7- 6e , r = 16e , (B18) 


with the slow roll parameters evaluated along the infla¬ 
tionary trajectory 




|VF|2 


2V 


r^ = M. 


2 ^^V V,J VjV 

V 


m2 ^4 4. (B19) 

^ 101 101 


The dotted line in Figure shows the power spectrum 
for TZ obtained through the slow-roll relation (|B18 ), with 
the slow roll parameter e evaluated through (B19i. This 


differs by about 2% from the exact result indicated by 
the solid line. The discrepancy is consistent with the 
accuracy of the slow-roll approximation (in this run, r] ~ 
1.5%). 


Appendix C: An incorrect simplification {tli = 0) 


In this Appendix we warn against a simplification that 
in general provides incorrect phenomenological results. 
The transformation Q rotates (up to negligible O 
corrections) the original axions 9 and (p of the starting 
potential § into the linear combinations (p and ip. Such 
fields are, respectively, the light and heavy directions of 
the inflaton potential close to the minimum at 9 = p = 0. 
For this reason, the inflationary evolution next to the 
origin occurs along the (p direction. One therefore may 
be tempted to simply set ip = 0 all throughout inflation. 
This would lead to the single field potential (11) with 
Ip = 0. 

In Figure[^we show with dashed lines the phenomeno¬ 
logical results obtained under this incorrect assumption. 
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r 



FIG. 10. Solid lines; correct phenomenological results; dashed 
lines: incorrect results obtained from incorrectly setting ip = 
0. The labels on the curves correspond to the following values 
of {rf, ta}: a = {1.5, 3} ; b = (1.5, 0.05} ; c = (1.65, 0.05} ; 
d = {1.8, 0.05} . For reference, the prediction of Natural In¬ 
flation is shown as the black dotted line. 


for four different choices of parameters. As a comparison, 


we also show with solid lines the correct phenomenolog¬ 
ical results, obtained from an exact numerical evolution. 
We see that the approximated results shown are in rea¬ 
sonable agreement with the correct ones only in the top 
part of the lines, where r is greatest. 

Each line is obtained by varying only and keep¬ 
ing the other parameters constant. The common point 
of all the lines corresponds to /^ —)■ oo. In this limit, 
^ only probes the immediate vicinity of the min¬ 

imum, where the model becomes identical to massive 
chaotic inflation (in practice, only the first quadratic 
term of the Taylor expansion of the potential around the 
minimum matters). As setting ■i/i = 0 is a good approx¬ 
imation close to the minimum, the approximate dashed 
lines approach the correct solid lines in that limit. How¬ 
ever, as decreases, the system probes the potential far 
from minima, where the ip = 0 assumption is incorrect. 
This explains why setting ijj = 0 does generally not pro¬ 
vide the correct phenomenological results for the model. 
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